Energy loss and dynamical evolution of quark px spectra 
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Average energy loss of light quarks has been calculated in a two stage equilibrium scenario where 
the quarks are executing Brownian motion in a gluonic heat bath. The evolution of the quark pr 
spectra is studied by solving Fokker-Planck equation in an expanding plasma. Results are finally 
compared with experimentally measured pion pr spectrum at RHIC. 
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Production of high transverse momentum [pt) parti- 
cles in heavy ion colhsion, in recent years, has assumed 
special interest. This is related with the phenomenon of 
high Pt particle suppression - dubbed as jet quenching. 
Actually energetic partons while passing through plasma 
lose energy which degrades the population of high pT 
hadrons. Experimentally such suppression has been ob- 
served at Relativistic Heavy Ion Collider (RHIC) at high 
Pt domain [1]. This phenomenon of 'jet quenching' can 
be used to extract the properties of early stage of the 
plasma, temporarily produced in the high energy heavy 
ion collision [2]. 

To study the modified pT spectrum of hadrons it is 
essential to estimate parton energy loss in the thermal 
bath of quarks and gluons. Partons in a plasma can dis- 
sipate energy in two ways, either by two body scatter- 
ing (collisional loss) or via the emission of gluons (radia- 
tive loss). Significant progress has been made in recent 
years to calculate partonic energy loss [2-9]. In many 
of these calculations, first, path length dependent energy 
loss is estimated, which, consequently is used to mod- 
ify the fragmentation functions that depopulate high pT 
hadrons [11]. Present approach is different in the sense 
that here we dynamically evolve the quark pT spectra 
for a given initial distribution. Importance of this has 
already been discussed in rcfs. [3,12] where the focus has 
been on the radiative loss. Even though this is the main 
mechanism of energy loss, under certain kinematic con- 
ditions, collisional loss could be comparable, or even be 
more, to its radiative counterpart [9]. This is particularly 
so in case of heavy quark because of the dead cone effect 
[10]. 

To estimate energy loss, we inject quarks with a given 
energy distribution and study the broadening of the same 
as the system expands and cools. In the present model 
gluons are thermalized at a much smaller time scale than 
quarks. Such two stage equilibration scenario, albeit in 
a different context, was considered sometime ago [13,14]. 
Possibility of earlier thermalization (even faster) of (soft) 
gluons is further accentuated, if the gluons are assumed 
to be in color glass condensate (CGC) state initially. In 



tion time scale for RHIC and LHC energies are estimated 
to be = 1.40 GeV-\ t, = 0.62 GeV"^ respectively [15]. 

Under such simplifying scenario, it thus reduces to the 
problem where quarks are executing Brownian motion 
in the gluonic heat bath. The evolution of the quarks 
are, therefore, governed by the Fokker-Planck equation 
(FPE). So we avoid solving full Boltzmann kinetic equa- 
tion (BKE) and approximate relevant collision integral 
in terms of appropriately defined drag and diffusion co- 
efficients [16-18]. 

The drag and diffusion co-efficients are calculated us- 
ing techniques of finite temperature perturbative quan- 
tum chromodynamics. The former is related to the quark 
energy loss, while the latter is related to the square of 
the momentum transfer [16,17]. Both of these quanti- 
ties show infrared divergences as the collisions are dom- 
inated by soft (t-channel) gluons. To cure this problem, 
we screen the interaction via hard thermal loop (HTL) 
corrected propagator which makes long range Coulomb 
interaction finite. 

To arrive at the relevant FP equation from BKE we 
assume that there is no external force and therefore. 



d_ 

dt 



TJT+Vp-Vr /(p,X,t)=C[/(p,X,i)] 



(1) 



Here, quarks have a phase space distribution which 
evolves in time and the collision term is evaluated by 
considering ultra-relativistic scattering of the quarks and 
gluons which eventually are expressed in terms of trans- 
port coefficients. Considering that the system expands 
in the longitudinal direction Eq.(l) takes the following 
form [19]: 



df{v,z,t) 
dt 



df {p,z,t) 
dz 



= C[f{p,z,t)]. 



(2) 



Here, Vpz — Pz/Ep (for light partons Ep — \p\). This 
equation can be simplified further for the central rapidity 
region which is boost invariant in rapidity, which implies 



/(PT,P2,z,i) = /(pt,P2,t). 



(3) 



the proper time. Using the Lorentz transformation rela- 
tion dT/dz\2=o = 0, 7^=0 = 1 and dp'^/dz\z=o = -p/t, 
one finds 



d.f Pz df 



dz 



t dpz 



(4) 



Therefore the Boltzmann equation takes the following 
form 

dfiPT,Pz,t). ( d Pz d \ 



§i ~ -^^P^'^^'*) = C'[/(pT,P.,t)]. (6) 

Evidently in Eq. 6, the second term on the left hand 
side represents the expansion while the right hand side 
characterizes the collisions. The latter can be written in 
terms of the differential collision rate Wp q 

C[/(pT,Pz,i)]= j d'q[Wj>+,,J(p + a,)-Wj>,J{p)] (7) 

which quantifies the rate of change of the quark momen- 
tum from p to p — q, VFp^q = (ir(p, q)/d^g, where F 

represent scattering rates. 

In a partonic plasma, small angle collisions, with para- 
metric dependence of 0{g^T), are more frequent than the 
large angle scattering rate. The latter goes as ^ 0{g'^T). 
Therefore the distribution function does not change much 
over the mean time between two soft scatterings. This 
allows us to approximate /(p + q) — ./(p)- In contrast, 
Wp.q, being sensitive to small momentum transfer, falls 
off very fast with increasing q. Therefore, we write 

W^P+q,q/(p + q) - W-P,q/(P) + 9i^(H^p.q/) 

With these approximation, Eq. 6 can be written as 



1 d 



[B,j{p)f{p)], (9) 



2 dpidpj 

where we have defined the following kernels, 

Ai = J d?qWp,^qi (10) 

Bij = j d^qWp^^qtqj (11) 

^ = ^ I f^?$<'«IA^lLo/W(i + /(*')) 



(fk'd^kcfq, .,_,.2 ,,,, 



(1 + f{k'))6'iq_ -~k' + k)5{u - Vk' • q)S{u; - Vk • q). 

First we consider Compton scattering {gq{q) — > gq{q)) 
for which 



9 us 

In terms of explicit momentum variables it reduces to, 
p'k'pk 



(14) 



\M\Uo 



(g2 - a;2)2 



{1 - cos9pk){l - cosOp'k')- (15) 



We label the incoming four momenta of the test and bath 
particles as P and K respectively and the corresponding 
final momenta are P' and K' . In Eq.(14) s,u and t are 
Mandelstam variables with t = uj'^ — q^ and s = —{u + t). 
The scattering angles, following ref. [20], can be written 
as 



^ ui t 

COSOpg =-+— 

q 2pq 

UJ t 

COsOkq = 



2pq 



(16) 



The parton parton collision is dominated by the soft 

momentum transfer. Therefore, we consider small an- 
gle scatterings to estimate the leading contribution, In 
this limit t ^ and s = —u. Hence the collision is 
dominated by the t channel, \M-\'^ g'^2s^ jt^ . We also 
have cos9kq ~ cosOk'q ^ ujjq ^ cosOpq ~ cosOpiq and 
cosOpk — uo'^/q^ -I- (1 — Lo'^ /q^)coscf)pk- With these az- 
imuthal angle averaged matrix element becomes, 



(|A4|)^o^ I2gyk'pk{l - u?lq^fl{q'^ - J^f 

~ 125VfcV9' (17) 

It might be noted that Eq. 12 is symmetric under 
k <-> A;', which allows one to write, 

. , 2^ V f cPk' d^k d^q , , , ..n 

l[f[k){i + f{k')) - f{k'){i + f{k))] 

<53(q - k' + k)(5(w - Vk' • q)<5(w - Vk • q). (18) 



recognizing the fact that lj is small compared to the mo- 
mentum, we have 



Furthermore, 



= j 2kf{k)dk = 



(20) 



With the help of these identities, the drag coefficient can 
easily be calculated at the leading log order: 



6p 



JC-, 



(21) 



where C = J ^ [15]. Evidently, C is infrared singular. 
Such divergences do not arise if close and distant col- 
lisions are treated separately. For very low momentum 
transfer the concept of individual collision breaks down 
and one has to take collective excitations of the plasma 
into account. Hence there should be a lower momen- 
tum cut off above which bare interactions might be con- 
sidered. While for soft collisions medium modified hard 
thermal loop corrected propagator should be used [8,9]. 
ft is evident that Eq. 12 actually gives dE/dt or the en- 
ergy loss rate [9] that can be related to the drag coeffi- 
cient. 

Bij can be decomposed into longitudinal and trans- 
verse components: 



R _ R a PiPi \ , R PiPj 



(22) 



Explicit calculation shows that the off diagonal compo- 
3ij vanish. 

i^g'^ f d^kd^qduj 



nents of Bij vanish 



Bt 



(2 



■l_ f 



2k2k'2p2p 



-6{lj - Vp • q)(5(w - Vk • q) 



(23) 



Here, in the small angle limit qi ^ lo and qt — \f<^ 
With all these, in the leading log approximation 



B, = '-T^jC 



It' 



(24) 



(25) 



Here the fluctuation-dissipation theorem, is found to be 
satisfied automatically. This connects drag and momen- 
tum diffusion constants giving rise to Einstein's relation 
B{p^) = 2TEA{p^), where B,j{p^) = 5ijB{p'^). To re- 
move arbitrariness related to the momentum cutoff scales 
hidden in £, hard thermal loop corrected propagator is 
used [9]. 

In the coulomb gauge, we can define Dqo = and 



/^l{qo,q)-^ =q^ -^Lol 



— In 2 

q qo-q 



(26) 



^t{qo,q) ^ =ql 



q +2'^p 



90 (go 



qo 



ql 



2g3 



qo 



With this, matrix element in the limit of small angle 
scattering, we get the following expression for the squared 
matrix element. 



\M\^ = 9^C,,l(!>{EE,f\I\i{q^,q) 
+ {vx q).{vi X q)At{qo,q)f 



(28) 



with V = p, vi = pi and Cqq is the color factor. With 
the screened interaction, the drag and diffusion constants 
can be calculated along the line of ref . [9] . 

Now to calculate average energy loss of the light quarks 
in an expanding partonic plasma we inject test quarks 
having following distribution at time ti, 



f{PT,Pz,t = ti 



Nd'^ipT 



■ Pto)S{Pz - Pzo) (29) 



Bjorken cooling law [21], TiT^^'^^ = tT^/^^ is used to de- 
scribe the space-time evolution of the plasma. Here Cg is 
the velocity of sound. 

As the time progresses the system expands and 
f{pT,Pz,t) evolves according to Eq. 9, which is solved 
numerically. The average energy loss is given by 



{AE) = (E) - Eo 



(30) 



where {E{t)) = J d^pEf{pT,Pz,t)/ J d^pf{pT,Pz,t). Eo 
corresponds to the peak of the initial spectrum, i.e. 
Eo =< E{ti) >. 
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In Fig. 1, we present results for the energy loss as a 
function of time for various initial energies with means 
indicated in the legend. It might be noted that at early 
times fractional energy loss is independent of the average 
incoming energy. 

Next we investigate time evolution of quark pt spec- 
tra. Assuming that the high pT partons materialize into 
hadrons outside the system and the hadronization pro- 
cess does not affect the shape of the pt distribution dras- 
tically. Under these circumstances the pion pr spectra 
may be taken to be proportional to quark pT spectra for 
a given time. Therefore, we take the initial pt distribu- 
tion of quark to be proportional to the pion pt spectrum 
as measured in p-p collision [23] : 



fiPT,Pz,t = U) = 



No dNjy = 0) 
cPprdy 
1 



Pt 

No 



(31) 



where, = 9.97 and po = 1.212 and iVo(^^o) is the 
normalization constant. The final spectrum is obtained 
by assuming the transition temperature, ~ 190 MeV 
at a time t ^ 2 fm/c. The results arc compared with 
the measured px spectrum of pions (for pT ^ 2 GeV) 
at RHIC [22] in Fig.2. The data is well described for 
Pt = 3 - 8 GeV. 



data 

t=0.3 fm/c 
t=lfin/c 




10" 



Ti = 350 MeV, ti = 0.3 fm/c 



3456789 10 

Pt (GeV) 

FIG. 2. Time evolution of the transverse momentum dis- 
tributions. 

In conclusion, in the present work, we estimate quark 
energy loss considering stochastic nature of the interac- 
tion by solving FP equation. Relevant drag and diffu- 
sion coefficients have been calculated in the soft collision 
limit. To highlight the importance of collisional energy 
loss radiative processes is excluded. We have analysed 
the pion pT spectrum measured at RHIC by the PHENIX 



collaboration [22]. The two body scattering is found to 
give reasonable amount of quenching required to explain 
the data. Recently transverse momentum spectrum for 
the D meson has also been measured via single electron 
Pt distribution [10,24]. The present formalism can be 
applied to analyse the data reported in ref. [24]. Such 
investigations are currently under progress [25]. 
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